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Abstract 

Project a collection of points on the high-dimensional sphere onto 
a random direction. If most of the points are sufficiently far from 
one another in an appropriate sense, the projection is locally close in 
distribution to the Poisson point process. 



1 Introduction 

Let xi,X2, ■ ■ ■ ,a;„ be n points on the {d — l)-dimensional sphere S'^~^. We 
assume that n, d, and the points themselves depend on an implicit parameter, 
so that d — i> oo, n — >■ oo. Consider the normalised projections (xj, \fdU) of 
the points onto a random direction U on the sphere. 

Fix a G M, and denote 

n 

= UU) = 5^ 5 ( • -n{{xj,^/dU) - a)) . 
i=i 

This is a point process, i.e. a random locally finite integer-valued Borel mea- 
sure on M. The (homogeneous) Poisson process with intensity A (A > 0) is a 
point process rj such that 

r/(5) ~ Pois(A -mesfi) (1) 
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for any Borel set B cM. The reader may find further properties of Poisson 
point processes in the book of Reiss [2|. 

Theorem. Assume that, for any e > 0, 

where the sum is over all pairs i < j such that \{xi,Xj)\ > e. Then, for any 
a ^'E., converges to the Poisson process with intensity 

0(a) = (27r)-i/=^exp(-aV2) , 
in the following sens^: for any hounded Borel set / (s 

UI) ^ Pois(0(a) 

That is, the random variables {xj,\fdU) behave locally as independent 
samples from the Gaussian distribution. 

Example. It is not hard to see that the condition ([2]) is fulfilled for the 

vertices of the discrete cube: 

n = 2'^ , Xj- = (±l/v^,--- ,±l/v^) . (3) 



2 Proof of Theorem 

The proof is based on the following elementary (and well-known) lemma, 
sometimes referred to as Archimedes' theorem. The lemma follows from the 
fact (discovered by Archimedes for d = 3) that the projection of the uniform 
measure on S'^~^ onto a ((i — 2)-dimensional subspace is the uniform measure 
on the unit ball of this subspace. 

Lemma 1. Let yi,--- ,yk be pairwise distinct points on the sphere S'^~^, 
k < d — 2. The projections Hj = {yj, \fdU) of yj on a random direction have 
joint density 

where Mjj> = {yj,yj')- 

^that is stronger than weak convergence in distribution 
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Let / (£ M, and fix A; G N. Denote by Nj the number of points of in I, 
= ^a{I)- Let us show that Nj converges in distribution to the Poisson 

law. 

Step 1: First, let us assume that 

max\{xj,Xj')\=o{l) . (4) 
37^3' 



Then proceed as follows: 
E 



^ ' l<ji<-<jfc<n^=l (5) 

= 5^Ej]l„+„-i,((x,,,v^f/)) . 

s 

Denote = {xj^, VdU), and let Mgs' = {xj^,Xj^,). According to Lemma [H 
the joint density of H = {Hi, ■ ■ ■ , H^) is equal to 

^^^^ " T{{d - k)/2) {ndf/^ "^""^ '^^^ '"^^^ ■ 

Now, according to (jl]), M = 1 + o(l), where the o(l) term tends to zero 
entry- wise and hence also in norm (recall that k is fixed). Thus 

T{d/2) 



where the o(l) term is uniform in h. Recalling that d 00 (whereas k is 
fixed), we see that 

'^'^= m%2) J)./. exp(-|/f/2)(l+o(l)) 
= (21)-"'- exp(-|A|V2) (1 + 0(1)) , 
uniformly on compact subsets of M.^. Therefore 

k 

E n la+n-i/(^s) = ilia + n-h)f (1 + 0(1)) , 
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where 7 = A^(0, 1) is the standard Gaussian measure. The set / is bounded 
and fixed, whereas n — 00, hence 



El[l,+^-ij{Hs) = Uia)"^) (l + o(l)). 

s=l ^ ^ 

Returning to ([5]), we deduce: 

That is, the factorial moments of Nj tend to those of the Poisson distribution 
Pois(0(a) mes/). The Poisson distribution has (better than) exponential 
tails, thus Nj converges in distribution to Pois(0(a) mesJ). 

Step 2: Now let us relax the assumption (jl]). First, ([2]) implies that one can 
choose e — ^ so that 

^'min I — = ^ n] = o(n^) . (7) 

Let 

pi{x, x') = P I {x, \fdU), (x', \fdU) G a + n'^I^ . 

Lemma 2. x') < Cj min — | (x, x') l)^"*^/^, n^-*^) . 

Proof. By Lemma [Tj the joint density of (x, y/dU), (x', \fdU) is given by 



d-4 
2 



where 

jVf =1 ^ 
\^ (x, x') 1 

Therefore 

det-i/'M = (1 - (x,xy)"'/' = 0(1) ■ (1 - |(x,x')|)"'/' 
andp(/i) = 0(1)(1- |(x,x')|)"^/^. Thus 

Pj(x,x')= !! p{h)dh<Cin-^{l-\{x,x')\)-^'^ 

J J{a+n-^I)^ 



Also, 

x') <F {x e a + n-^l} < Cin-^ . 

□ 

According to the lemma and ([7]), W{A) = o(l), where 

A = |3j 7^ / I \{xj,Xj>)\ > e, {xj, VdU), {xy, VdU) e a + . 

Repeating the argument of Step 1, we see that the conditional distribution of 
Nj given -lA (the negation of A) tends to Pois(0(a) mes/). Thus the same 
is true for Nj itself. 

□ 

3 Some remarks 

1. Diaconis and Freedman [1] have proved the following: if, for any e > 0, 

#{j|||x,f -l|>£}=o(n), (8) 
i^{j,k\\{xj,Xk)\> e} =o{n'^) , (9) 

then the empirical distribution 

n-^ J26{»-{xj,VdU)) 

j 

converges (weakly, in distribution) to the standard Gaussian law. Our result 
can be seen as a local version of this statement. 

2. The conditions ([8])- ([9]) are not sufficient for the conclusion of our theorem, 
as one can see from the following example: 

{xi,X2,--- ,Xn) = (ei,--- ,ed,ei,--- ,eisdi) 

(where (ei, ■ ■ ■ , e^) is the standard basis in M°', and < 5 < 1 is an arbitrary 
constant.) 

3. The assumption that Xj G S"^^^ in our theorem can be relaxed. For 
example, the Diaconis-Freedman assumption ([H]) is sufficient for a 7^ 0. 
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4. For any S > 0, one can construct a S-net on S'^~^ for which the assumption 
([2D is satisfied. Indeed, if the distribution of the points in the net is sufficiently 
regular, 

^ ^n^ [[ dad{x)dadiy) ^ ^^^2^ 

V^l - \ {Xi,Xj)\ Jj\{x,y)\>e \/l - \ {x,y)\ 

as d ^ 00. 

5. The theorem and the proof can be easily extended to random projections 
onto an r dimensional subspace, where r is any number (fixed, or slowly 
growing with d). 

6. One may ask whether it is possible to reduce the randomness in the 
conclusion of the theorem, and still have (at least, weak) convergence to the 
Poisson process. For example, one may project the point xj onto a random 
Bernoulli direction B = (±1, ■ ■ ■ ,±1). Even for the points ([3]), the limit 
will not be Poisson, since all the projections will be integer multiples of 
1/Vd. Instead, one can consider a random perturbed Bernoulli direction: 
Bs = (±(1 + ei),--- ,±(1 + £d))- Is it true that the projections of the 
points ([3]) are asymptotically Poisson for a 'generic' perturbation e? Is there 
a natural arithmetic condition on e that ensures that the projections are 
asymptotically Poisson? 

7. It may also be interesting to consider projections of points {xj} for which 
the condition ([2]) is violated. Which point processes can appear in the limit, 

as n,d ^ 00? 
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